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ABSTRACT 

Context. Wide-field surveys are a commonly-used method for studying thousands of objects simultaneously, to investigate, e.g., the 
joint evolution of star-forming galaxies and active galactic nuclei. Very long baseline interferometry (VLBI) observations can yield 
valuable input to such studies because they are able to identify AGN unambiguously in the moderate to high-redshift Universe. 
However, VLBI observations of large swaths of the sky are impractical using standard methods, because the fields of view of VLBI 
observations are of the order of 10" or less, and have therefore so far played only a minor role in galaxy evolution studies. 
Aims. We have embarked on a project to carry out Very Long Baseline Array (VLBA) observations of all 96 known radio sources in 
one of the best-studied areas in the sky, the Chandra Deep Field South (CDFS). The challenge was to develop methods which could 
significantly reduce the amount of observing (and post-processing) time, making such a project feasible. 

Methods. We have developed an extension to the DiFX software correlator which allows one to efficiently correlate up to hundreds of 
positions within the primary beams of the interferometer antennas. This extension enabled us to target many sources simultaneously, 
at full resolution and high sensitivity, using only a small amount of observing time. The combination of wide fields-of-view and high 
sensitivity across the field in this survey is unprecedented. 

Results. We have observed with the VLBA a single pointing containing the Chandra Deep Field South, in which 96 radio sources were 
known from previous observations with the Australia Telescope Compact Array (ATCA). From our input sample of 96 sources, 20 
were detected with the VLBA, and one more source was tentatively detected. The majority of objects have flux densities in agreement 
with arcsec-scale observations, implying that their radio emission comes from very small regions. Two objects are visibly resolved. 
One VLBI-detected object had earlier been classified as a star-forming galaxy. Comparing the VLBI detections to sources found in 
sensitive, co-located X-ray observations we find that X-ray detections are not a good indicator for VLBI detections. 
Conclusions. We have successfully demonstrated a new extension to the DiFX software correlator, allowing one to observe hundreds 
of fields of view simultaneously. In a sensitive observation of the CDFS we detect 21 % of the sources and were able to re-classify 7 
sources as AGN which had not been identified as such before. Wide-field VLBI survey science is now coming of age. 

Key words. Techniques: interferometric, Galaxies: active, Galaxies: evolution 



1. Introduction 

VLBI observations provide the highest resolution in observa- 
tional astronomy, but at the cost of tiny fields of view. Until 
recently, only objects with brightness temperatures of the order 
of 10 6 K and greater could be observed, and since such objects 
are distributed sparsely on the sky, the small field of view has 
not been much of a limitation in the past. However, recent ad- 
vances in technology allow the use of much larger bandwidths 
in the observations, allowing much fainter objects to be detected. 
Their density on the sky is so high that many can be found in the 
antenna primary beams no matter where one looks. 

Unfortunately, the long baselines of VLBI observations tra- 
ditionally allow only objects within a few arcseconds of the 
phase centre to be observed, before effects arising from time 
and frequency averaging wash out the signal. The bottleneck has 
been correlator capacity (and post-processing power). Software 
correlators, and in particular DiFX (Deller et al. 2007 1, are much 
more flexible than traditional correlators, giving observers the 



possibility to image objects anywhere in the antenna primary 
beams. We have developed the necessary software, expertise and 
techniques to image all known radio sources within the primary 
beam of the interferometer elements using a single VLBI obser- 
vation, the results of which we present in this paper. 

There are numerous scientific motivations for carrying out 
wide-field VLBI surveys. One of the most pressing is to identify 

(1998) 



active galactic nuclei (AGN). For example, IHeisler et al. 



obtained VLBI detection rates of % for starbursts compared to 
more than 90 % for AGNs in their sample of infrared-selected 
galaxies. However, as instrumental sensitivities increase, VLBI 
observations have become more sensitive to detection of ra- 
dio supernovae in the relatively nearby Universe, such as those 
|1998b| l, and so a VLBI detection no 



in Arp220 ( Smith et al 



longer unambiguously implies an AGN. For example, |Kewley| 



et al. ( 2000[ ) observed a sample of 61 luminous infrared galax- 
ies and found compact radio cores in 37 % of starburst galaxies 
and 80% of AGNs. However, the radio luminosity distribution 
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of the Kewley et al. compact cores showed a bimodal distribu- 
tion. Most higher radio luminosity cores (> 2 x 10 21 WHz -1 ) 
are AGNs, while most lower radio luminosity cores (< 2 x 
10 21 WHz [ ) are starbursts. We conclude that high luminosity 
radio cores (» 2 x 10 21 WHz ') VLBI detections are almost 
certainly AGN, whereas lower luminosity cores may be caused 
either by AGN or by supernova activity (also |Smith et al.|19 98a). 
For our detection limit of ~ 0.5 mJy, we can be confident that 
VLBI detections above z > 0.1 represent AGNs. Below this red- 
shift, a VLBI detection is an ambiguous indicator of AGN ac- 
tivity. Furthermore, the amount of radio emission from an AGN 
does not correlate well with emission at other wavelengths; in 
particular there appears to be no correlation between a galaxy's 
nuclear radio and its large-scale FIR emission (Corbett et al. 
[20021 . 

It is now recognised that AGN play an important role in star 
formation and galaxy evolution. They are not only found in pow- 
erful quasars and radio galaxies, but also in the loc al universe, 



Filho et al. 



with radio luminosities as small as 10 20 W Hz 1 (e.g. 
2006). AGN feedback can either push back and heat the gas, re 
ducing the formation of stars in the galaxy, or compress the gas 
clouds, triggering the same process. Particularly at high redshift, 
both AGN and star formation processes are likely to be impor- 
tant in a large fraction of galaxies ( |Bower et al. |2006| l, but we 
can not yet measure the fraction of the luminosity (bolometric 
and radio) generated by each process, nor how they are influ- 
enced by feedback. Whilst observations at other wavelengths, in 
particular in the optical and X-ray regime, can identify accretion 
onto black holes, these data cannot reliably indicate whether or 
not an AGN also produces radio jets which interact with their 
surroundings. 

Furthermore, VLBI observations of a substantial number 
of galaxies yield information about specific classes of object. 
For example, the compact steep-spectrum (CSS) and gigahertz- 
peaked spectrum (GPS) sources make up a significant fraction 
(10% to 30 %, |Q'Dea]|T99"8l of radio sources. CSS and GPS 
sources are very small, yet strong radio sources with either 
steep spectra (a < -0.5, S oc y") or spectra that peak be- 
tween 500 MHz and 10 GHz. Since only bright members of this 
class have yet been studied with VLBI observations, a wide-field 
VLBI survey can yield data in a previously inaccessible flux den- 
sity regime. Another example is the infrared-faint radio sources 
(IFRS), a mysterious class of object characterised by strong (up 
to 20mJy) 1.4 GHz radio emission but a flux density of less than 
~ 1 /iJy in the near-IR at 3.6 /im (Norr is et al.||2006|l. It is now 
established that they host AGN ( |Norris et al.|2010| |Middelberg 
|et al.|[20T0| and references therein), but their relation to other 
classes of object is unclear. 

The connection between radio and X-ray emission has re- 
cently been investigated by Dunn et al. ([2010), using a sam- 
ple of nearby (D < 100 Mpc) X-ray bright elliptical galaxies. 
They find that almost all sources are also radio emitters and 
most display a likely interaction between the X-ray and radio 
emitting plasmas. On the other hand, they find no correlation be- 
tween radio and X-ray luminosity. For a study of AGN which are 
not recognized as such in optical observations ("elusive AGN"), 



Maiolino et al. ( 2003 \ used a VLBI-selected sample of 18 galax- 



ies to maximise the probability of X-ray detections. They find 
that these objects are heavily obscured (column densities ex- 



ceeding 10 cm ) and that their space density is comparable to, 
or exceeding that of, optically classified Seyferts. We note how- 
ever, that both works are difficult to compare to what we present 
here. Dunn et al. (2010) used VLA observations of nearby el- 
liptical galaxies, and even though the sample by Mai olino et al.| 



(2003) was VLBI-selected, the parent samples were chosen for 
high IR luminosity and detectability with VLBI. In contrast, our 
sample is only flux density-limited, and we make no further se- 
lection. 

Throughout the paper, we adopt a flat ACDM cosmology 
with H = 71.0, Qm — 0.27, Q vac = 0.73. 

2. The target sources 

Executing a pilot widefield VLBI observation such as the one 
described below required the development of new processing 
techniques, their implementation in software and the purchase of 
necessary computing infrastructure. A perfect test-bed for these 
techniques was the Chandra Deep Field South (CDFS). It has 
rich complementary coverage at many wavelengths, particularly 
in the radio regime. However, even with novel processing tech- 
niques it was not possible to image the entire area contained 
within the antennas' primary beams, nor would this be a promis- 
ing exercise, since any flux detected in the VLBI observations 
must have been detected in existing compact-array interferome- 
try data because of their sensitivity to lower-brightness temper- 
ature sou rce j*| We the refore used the ATCA observations pub- 
lished by Norris et al. (2006) as an input catalogue, and aimed 
at imaging all 96 targets from that catalogue in a region cen- 
tred on the GOODS/CDFS and contained within the VLBA pri- 
mary beam area. An overview of the observed field is given in 
Figure[T[ and the details of the target sources are listed in Table[T] 

3. Observations and data analysis 

3.1. Observations 

We observed on 3 July 2007 the GOODS/CDFS area with 
the VLBA at 1383 MHz, centred on RA 03:32:34.0392, 
Dec -27:50:50.748 (J2000). Two polarisations were recorded 
across a total bandwidth of 64 MHz, using 2-bit sampling, re- 
sulting in a recording bit rate of 512 Mbps. The fringe finders 
3C 454.3 and 4C 39.25 were observed in intervals of 2.5 h, and 
the phase calibrator NVSS J034838-274914 was observed for 
1 min after each 5 min scan on the target field. The elapsed time 
of the observations was 9 h, but the low declination of the target 
field and the calibrator observations reduced the total integra- 
tion time on the target field to 178 baseline-hours. According 
to the VLBA observational status summarj|^] the baseline sen- 
sitivity of the VLBA at 1.4GHz is 3.3 mJy in a 2 min observa- 
tion using 256 Mbps recording, or 0.43 mJy per baseline-hour 
using 512 Mbps. This sensitivity results in a thermal noise limit 
of our images of 32 yuJy beam -1 . Accounting for the loss of about 

40 % of the data in self-calibration on our in-beam calibrator 
S503 (see Table [TJ, the estimate has to be increased to about 

41 /iJy beam -1 . However, these estimates are typical for observa- 
tions near the zenith, but our observations were carried out pre- 
dominantly at low elevations (at an average elevation of around 
20°), where the system temperatures were significantly higher, 
resulting in a lower expected system sensitivity. For example, 
the average r sys of the Fort Davis station when observing the 
target field was 27 % higher than compared to a fringe finder ob- 
servation at 75° elevation. At Kitt Peak, the increase was 32 %, 
and at Los Alamos 42 %. The final on-axis noise was found to be 
55 yuJy beam -1 , only 27 % higher than at an average declination 
of 75°. 



This reasoning ignores variability and transient events, see Lenc 



et al. 1 2008 1 for an example 
/ http://www.vlba.nrao.edu/astro/obstatus/current/obssum.html 
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Table 1: The 96 sources in our sample taken from |Norris et al. (2006). Columns are the ID we use in this paper, the IAU designation, 
right ascension and declination, and ATCA flux density in mJy. 



ID 


IAU desij 


;nation 








RA (J2000) 


Dec (J2000) 


5 1 .4GHz/mJy 


S329 


ATCDFS 


J033123 


30- 


•274905 


6 


03 


31 


23 


305 


-27 


49 


05 


630 


1 


1 


S331 


ATCDFS 


J033124 


89- 


-275208 


3 


03 


31 


24 


892 


-27 


52 


08 


320 


35 


5 


S339 


ATCDFS 


J033127 


51- 


•274440 


1 


03 


31 


27 


519 


-27 


44 


40 


130 





1 


S340 


ATCDFS 


J033128 


59- 


•274934 


9 


03 


31 


28 


594 


-27 


49 


34 


930 





5 


S343 


ATCDFS 


J033130 


07- 


•275602 


5 


03 


31 


30 


077 


-27 


56 


02 


520 





9 


S347 


ATCDFS 


J033130 


72- 


-275734 





03 


31 


30 


720 


-27 


57 


34 


050 





3 


S359 


ATCDFS 


J033138 


47- 


-275942 





03 


31 


38 


473 


-27 


59 


42 


060 





4 


S360 


ATCDFS 


J033139 


57- 


-274119 


4 


03 


31 


39 


572 


-27 


41 


19 


400 





2 


S367 


ATCDFS 


J033146 


08- 


-280027 


2 


03 


31 


46 


080 


-28 


00 


27 


260 





2 


S368 


ATCDFS 


J033146 


57- 


-275734 


8 


03 


31 


46 


572 


-27 


57 


34 


840 





3 


S374 


ATCDFS 


J033149 


87- 


-274838 


8 


03 


31 


49 


876 


-27 


48 


38 


800 


2 


1 


S376 


ATCDFS 


J033150 


07- 


-273947 


1 


03 


31 


50 


076 


-27 


39 


47 


190 





6 


S377 


ATCDFS 


J033150 


78- 


-274704 





03 


31 


50 


788 


-27 


47 


04 


030 





5 


S380 


ATCDFS 


J033152 


12- 


-273926 


4 


03 


31 


52 


122 


-27 


39 


26 


450 





8 


S382 


ATCDFS 


J033153 


41- 


-280221 


1 


03 


31 


53 


412 


-28 


02 


21 


130 





5 


S393 


ATCDFS 


J033201 


44- 


-274647 


5 


03 


32 


01 


440 


-27 


46 


47 


590 


49 


1 


S395 


ATCDFS 


J033202 


83- 


-275612 


6 


03 


32 


02 


836 


-27 


56 


12 


610 





2 


S396 


ATCDFS 


J033203 


84- 


-275804 


6 


03 


32 


03 


841 


-27 


58 


04 


610 





2 


S403 


ATCDFS 


J033208 


54- 


-274647 


7 


03 


32 


08 


548 


-27 


46 


47 


750 





2 


S404 


ATCDFS 


J033208 


67- 


-274734 


4 


03 


32 


08 


670 


-27 


47 


34 


430 


1 


7 


S405 


ATCDFS 


J033209 


75- 


-274247 


4 


03 


32 


09 


753 


-27 


42 


47 


420 





2 


S407 


ATCDFS 


J033210 


13- 


-275936 


8 


03 


32 


10 


132 


-27 


59 


36 


830 





9 


S410 


ATCDFS 


J033210 


79- 


-274627 


8 


03 


32 


10 


793 


-27 


46 


27 


850 





2 


S411 


ATCDFS 


J033210 


91- 


-274415 


2 


03 


32 


10 


915 


-27 


44 


15 


210 


2 


7 


S412 


ATCDFS 


J033211 


00- 


-274053 


8 


03 


32 


11 


002 


-27 


40 


53 


860 





3 


S414 


ATCDFS 


J033211 


64- 


-273726 


1 


03 


32 


11 


645 


-27 


37 


26 


160 


3 


6 


S415 


ATCDFS 


J033213 


07- 


-274351 





03 


32 


13 


077 


-27 


43 


51 


070 


1 


2 


S416 


ATCDFS 


J033213 


27- 


-274241 


3 


03 


32 


13 


277 


-27 


42 


41 


300 





1 


S417 


ATCDFS 


J033214 


17- 


-274910 


7 


03 


32 


14 


179 


-27 


49 


10 


780 





1 


S418 


ATCDFS 


J033214 


89- 


-275640 


7 


03 


32 


14 


895 


-27 


56 


40 


790 





2 


S421 


ATCDFS 


J033217 


05- 


-275846 


5 


03 


32 


17 


051 


-27 


58 


46 


570 


2 


6 


S423 


ATCDFS 


J033218 


01- 


-274718 


5 


03 


32 


18 


011 


-27 


47 


18 


570 





4 


S425 


ATCDFS 


J033219 


17- 


-275407 





03 


32 


19 


177 


-27 


54 


07 


040 


8 


8 


S427 


ATCDFS 


J033219 


46- 


-275219 


1 


03 


32 


19 


465 


-27 


52 


19 


120 





3 


S429 


ATCDFS 


J033221 


07- 


-273529 


7 


03 


32 


21 


074 


-27 


35 


29 


720 





2 


S430 


ATCDFS 


J033221 


22- 


-274435 


3 


03 


32 


21 


223 


-27 


44 


35 


310 





2 


S432 


ATCDFS 


J033222 


59- 


-280023 


6 


03 


32 


22 


590 


-28 


00 


23 


690 





4 


S433 


ATCDFS 


J033222 


62- 


-274426 


5 


03 


32 


22 


627 


-27 


44 


26 


550 





1 


S434 


ATCDFS 


J033223 


70- 


-273649 


1 


03 


32 


23 


701 


-27 


36 


49 


180 





1 


S435 


ATCDFS 


J033223 


79- 


-275844 


9 


03 


32 


23 


793 


-27 


58 


44 


960 





2 


S436 


ATCDFS 


J033226 


82- 


-280453 


1 


03 


32 


26 


828 


-28 


04 


53 


190 





2 


S437 


ATCDFS 


J033226 


97- 


-274106 


7 


03 


32 


26 


975 


-27 


41 


06 


710 


16 


6 


S439 


ATCDFS 


J033228 


76- 


-274619 


7 


03 


32 


28 


763 


-27 


46 


19 


750 





3 


S440 


ATCDFS 


J033228 


79- 


-274356 


1 


03 


32 


28 


790 


-27 


43 


56 


100 


3 


8 


S442 


ATCDFS 


J033229 


84- 


-274423 


8 


03 


32 


29 


847 


-27 


44 


23 


880 





9 


S443 


ATCDFS 


J033229 


97- 


-274405 


4 


03 


32 


29 


974 


-27 


44 


05 


480 





4 


S444 


ATCDFS 


J033230 


52- 


-275911 


6 


03 


32 


30 


529 


-27 


59 


11 


680 





3 


S446 


ATCDFS 


J033231 


54- 


-280433 


5 


03 


32 


31 


540 


-28 


04 


33 


530 





3 



3.2. Correlation 

3.2.1. Bandwidth and time averaging effects 

Since the targets are scattered throughout the primary beam it 
was necessary to develop new correlation strategies to overcome 
the effects of bandwidth and time averaging. These effects are 
commonly called bandwidth smearing and time smearing, since 
they smear out emission from the target into the map and re- 
duce the observed amplitudes of sources away from the phase 
centre. Bandwidth smearing is regarded as the aperture synthe- 
sis equivalent to chromatic aberration, whereas time smearing 
can be regarded as similar to motion blur in a photograph when 
the shutter time is too long. An estimate of the magnitude of 
these effects can be found in Thompson et al. (2001). In our 



case, with sources having separations to the phase centre of up to 
15 arcmin, and long baselines of typically 5000km, amplitudes 
would be reduced to a few percent of their true values if a normal 
correlator setup with 500 kHz channels and 2 s integrations was 
used. 

To keep the amplitudes within 5 % of their true values, one 
would have to use a channel width of 4 kHz and 50 ms integra- 
tions. Such a correlation would result in around 3 TB of visibil- 
ity data, which would be very difficult to manage even on large 
general-purpose computers. Furthermore, current hardware cor- 
relators are not able to produce data with such high resolution. 
Therefore the only option was to use a software correlator, in 



particular the DiFX correlator developed by Deller et al. (2007 ). 
We used a development version of DiFX which included a new 



3 



Enno 



Middelberg et al.: Wide-field VLBA 
Table 1: 



Observations of 
(continued) 



the Chandra Deep Field South 



ID 


IAU designation 


RA 


Dec 


Si.4GHz/mJy 


S447 


ATCDFS 


J033232 


04- 


-280310 


2 


03 


32 


32 


041 


-28 


03 


10 


290 


23 


9 


S450 


ATCDFS 


J033233 


48- 


•275227 


2 


03 


32 


33 


480 


-27 


52 


27 


240 





1 


S453 


ATCDFS 


J033237 


73- 


-275000 


3 


03 


32 


37 


737 


-27 


50 


00 


330 





1 


S457 


ATCDFS 


J033238 


95- 


-275700 


6 


03 


32 


38 


954 


-27 


57 


00 


660 





1 


S458 


ATCDFS 


J033239 


57- 


-280312 


5 


03 


32 


39 


570 


-28 


03 


12 


570 





2 


S459 


ATCDFS 


J033241 


60- 


-280128 





03 


32 


41 


604 


-28 


01 


28 


060 





4 


S461 


ATCDFS 


J033241 


80- 


-280552 


5 


03 


32 


41 


808 


-28 


05 


52 


520 





3 


S462 


ATCDFS 


J033242 


83- 


-273817 


6 


03 


32 


42 


830 


-27 


38 


17 


660 


72 


3 


S463 


ATCDFS 


J033244 


16- 


-275142 


5 


03 


32 


44 


165 


-27 


51 


42 


580 





5 


S465 


ATCDFS 


J033245 


36- 


-280449 


7 


03 


32 


45 


363 


-28 


04 


49 


730 


2 


9 


S469 


ATCDFS 


J033245 


95- 


-275745 


2 


03 


32 


45 


956 


-27 


57 


45 


290 





2 


S472 


ATCDFS 


J033249 


19- 


-274050 


7 


03 


32 


49 


195 


-27 


40 


50 


710 


4 


7 


S473 


ATCDFS 


J033249 


33- 


-275845 


1 


03 


32 


49 


334 


-27 


58 


45 


100 





2 


S474 


ATCDFS 


J033249 


43- 


-274235 


3 


03 


32 


49 


430 


-27 


42 


35 


350 


3 


2 


S479 


ATCDFS 


J033252 


48- 


-275942 


1 


03 


32 


52 


480 


-27 


59 


42 


180 





2 


S481 


ATCDFS 


J033253 


34- 


-280200 


1 


03 


32 


53 


341 


-28 


02 


00 


150 





9 


S482 


ATCDFS 


J033256 


47- 


-275848 


4 


03 


32 


56 


472 


-27 


58 


48 


460 


1 


2 


S484 


ATCDFS 


J033257 


17- 


-280210 





03 


32 


57 


170 


-28 


02 


10 


030 


23 


1 


S485 


ATCDFS 


J033259 


22- 


-274325 


6 


03 


32 


59 


224 


-27 


43 


25 


640 





2 


S488 


ATCDFS 


J033302 


69- 


-275642 


4 


03 


33 


02 


692 


-27 


56 


42 


400 





2 


S489 


ATCDFS 


J033303 


29- 


-275326 


4 


03 


33 


03 


293 


-27 


53 


26 


470 





1 


S491 


ATCDFS 


J033304 


38- 


-273804 


3 


03 


33 


04 


381 


-27 


38 


04 


310 





1 
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multiple-field capability, and which was installed on a computer 
cluster at the Max-Planck-Institut fur Radioastronomie in Bonn. 



3.2.2. An extension to the DiFX correlator 



The DiFX software correlator is widely-used and verified (Deller 
et al. 2007 ). It replaced the hardware correlator of the Australian 
Long Baseline Array in 2008. Subsequent to our observations 
DiFX has also been adopted as the full-time correlator for the 
VLBA, and the multiple-field extension capability which is dis- 
cussed below has recently been incorporated into the VLBA im- 
plementation of DiFX and made available to the VLBA commu- 
nity. 

The implementation of "simultaneous multiple field centre 
correlation" is described in detail by |Deller et al.| ( [2~010[ ), and is 



briefly summarised here. The correlation is initially performed 
with high frequency resolution which is sufficient to minimise 
bandwidth smearing. Periodically, but still frequently enough to 
minimise time smearing, the phase centre of the correlated data 
is shifted from its initial location (which is usually the pointing 
centre) to a target source location. This shift requires rotating 
the visibility phases of each baseline by an amount equal to the 
difference in the geometric delay between the final and initial 
source directions, multiplied by the sky frequency. In effect, this 
corrects for the "unapplied" differential fringe rotation between 
the final and initial source directions. This phase shift is repeated 
for each desired source direction. After the phase shift is applied, 
the visibilities are averaged in frequency and continue to be av- 
eraged in time. Eventually, this results in an array of "normal- 
sized" visibility datasets, with one dataset per target source. The 
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3"33 m 30 s 33 m 00 s 33 m 30 s 33 m 00 s 

Right Ascension (J3000) 

Fig. 1: An overview of the observed area. The background image 
is a radio image of the CDFS made with the Australia Telescope 
Compact Array ( |Norris et al.||2006) . The rms of the image is 
around 20 /iJy beam 1 , and the faintest sources have flux densi- 
ties of lOOyuJy. The square indicates the ECDFS area observed 
with the Chandra satellite by |Lehmer et al.| ( [2005) > with an inte- 
gration time of 240 ks; the large circle indicates a typical VLBA 
antenna's primary beam size at the half -power level at 1.4 GHz; 
and the medium circle is the region covered uniformly with a 
2 Ms exposure with Chandra by Luo et al. (2008 1 (centred on 



the average aim point with a radius of 7.5', see their Figure 2). 
Crosses are drawn at the locations of the 96 targets taken from 
([2006 ) and small circles those targets which were 



Norris et al. 



detected with the VLBA 



3.3. Calibration 
3.3.1. Standard steps 

The calibration followed standard procedures used in phase- 
referenced VLBI observations, using the Astronomical Image 
Processing System, AIPS. Amplitude calibration was carried out 
using T sys measurements and known gain curves. Fringe-fitting 
was carried out on all calibrators to determine delays and phase 
corrections. The data of the phase calibrator were then exported 
to Difmap ( Shepherd|1997 1 for imaging. We have written an ex- 
tension to Difmap to allow export of the complex gains deter- 
mined in Difmap to AIPS for further processing and to apply 
the gains to the target data sets. Images were then made of the 
brighter sources to look for a source suitable for in-beam cali- 
bration. 

One of the targets, S503, was found to be bright enough 
for self-calibration. We applied the phase calibrator gains to the 
S503 data set and used two cycles of phase-only self-calibration 
to refine the gains, when convergence was reached. Amplitude 
self-calibration was found to converge on less than 30 % of the 
data and resulted in large gain fluctuations, and was therefore 
not used on S503. Furthermore, S503 was found to be resolved 
on long baselines, resulting in the number of useful visibilities 
reduced by » 40 % after self-calibration. The noise in images 
(before primary beam correction) was found to be of the order of 
55;uJy bearrT 1 . This sensitivity is in reasonable agreement with 
our expectations, given that (i) the number of visibilities was 
much reduced in self-calibration, and (ii) the system tempera- 
tures in our data were significantly increased compared to ob- 
servations near the zenith, an effect arising from the average an- 



tenna elevations being of the order of only 20 ° (see Section 3.1 

Since the phase response of a VLBA antenna is constant 
within the primary beam, and since the geometric delays had 
been taken care of at the correlation stage, the amplitude, phase, 
and delay correction could simply be copied from one target 
source to another, and an image could be formed. However, an- 
other correction is needed to compensate for the amplitude re- 
sponse arising from primary beam attenuation, which reduces 
the apparent flux density of a source by up to 50 %. 



field of view of each of these datasets is of the order of 13", 
at which point bandwidth and time smearing would reduce the 
observed amplitudes by 5 %. 

In comparison to the normal correlation operation, which re- 
quires operations on every baseband sample, at a time scale of 
~ns, these shift/average operations need only be applied rela- 
tively infrequently (10 ms is sufficient to shift past the edge of 
the primary beam with minimal time smearing), and thus adding 
many phase centres to the correlation is relatively computation- 
ally cheap. [Deller et al.] f2010) show that hundreds of phase cen- 
tres can be added within the primary beam at a cost only ~3 
times that of a single traditional correlator pass. Hence each in- 
dividual phase centre requires far fewer operations than a tradi- 
tional "multi-pass" correlation where the full correlation opera- 
tion is repeated for each phase centre. 

|Deller et"aT] ( |2010] l also demonstrate that the accuracy of the 
data obtained in this manner is on par with the accuracy of data 
obtained in the traditional way. Using a small test data set ob- 
served at a frequency of 8.4 GHz, they find the amplitude and 
phase differences of the two methods to be 0.09% and 0.014° 
on average, which is in agreement with effects expected to arise 
from numerical errors. 



3.3.2. Primary beam correction 

Unlike in compact-array interferometry, where the primary beam 
correction is carried out in the image plane, we have corrected 
for primary beam attenuation by calculating visibility gains. This 
is possible only because the fields of view are very small in our 
observations and so the attenuation due to the primary beam does 
not vary significantly across the image. In general the measured 
amplitude, A', of a visibility is related to the intrinsic amplitude, 
A, via 

A' - A -\Jgi(t, v)g 2 (t, v), (1) 

where g\{t,v) and g 2 (t, v) are the time and frequency- 
dependent gains of the two antennas of the baseline in the direc- 
tion of the target. The gains consist of many contributions from 
various instrumental and environmental effects, but the attenu- 
ation of the signal at each antenna arising from primary beam 
attenuation, b\{t, v) and b 2 (t, v), can be separated, yielding 

A' =A^b l (t,v)g' l (t,v)b 2 (t,v)g' 2 (t,v). (2) 

3 http://www.atnf.csiro.au/people/Emil.Lenc/tools/Tools/ 
Cordump.html 
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The gains g[(t, v) and g'Jt, v) are determined using r sys mea- 
surements and antenna gain curves, but the primary beam atten- 
uation is ignored in standard VLBI calibration procedures be- 
cause one usually observes a tiny field of view in direction of 
the optical axis where these terms can be ignored. 

Two effects need to be taken care of to determine b\(t, v) 
and btit, v) at any one time. First, in the case of the VLBA, 
the receivers are offset from the optical axis, and the beams 
of the right-hand and left-hand circular polarisations (RCP and 
LCP) therefore are offset on the sky. This effect is known as 
beam squint, and it scales with frequency (see |Uson & Cotton| 
2008 for a description of beam squint at the VLA). Second, the 



SC beam response 



correction factor = l/sqrt{response) 



beamwidths also are functions of frequency. Since the antennas 
are alt-azimuth antennas the sources rotate through the two off- 
set beams, and in general sources away from the pointing centre 
have different parallactic angles, because the separation of the 
antennas is large (i.e., at any given time antennas see the targets 
through different portions of their primary beams). This situation 
is illustrated in Figure [2] The angular separation (RCP-LCP) be- 
tween the two beams on the sky, averaged across all 10 antennas 
of the array, is -1.42' in azimuth and -0.58' in elevation (R. C. 
Walker, priv. comm.), or a total of 1.53', corresponding to 5.4 % 
of the average FWHM. This squint is in good agreement with 
theoretical predictions (Fom alont & Perley]|1999 1, and is small 
when the targets are close to the pointing centre (which is half- 
way between the two beams). However, near the FWHM point 
of an antenna beam the slope of the beam response is quite steep, 
and the difference between the RCP and LCP response can be up 
to 10%. 

The beamwidths and the amount and direction of beam 
squint have been measured at 1438 MHz (Walker, priv. comm.). 
We have scaled these values to the centre frequency of each IF to 
account for their frequency dependence. The relative locations of 
the antenna pointing centre, the target coordinates and the RCP 
and LCP beams on the sky were then computed in 1 min intervals 
for each IF, and the appropriate correction factor was calculated 
and saved to a new calibration table. These steps were carried 
out for each target data set separately. 

Primary beam corrections in VLBA observations are uncom- 
mon, and this is the first attempt to do so on a larger scale, as far 
as we are aware. Even though we have taken care to account for 
primary beam attenuation as accurately as possible, an ampli- 
tude error will unavoidably remain after correction. For exam- 
ple, the quadrupod legs and the subreflector blockage produce 
errors which have not been accounted for, and the primary beam 
shape and beam squint are likely to vary with elevation. We es- 
timate that this residual error is of the order of 10 %, which we 
add in quadrature to the other sources of error. 

3.3.3. Imaging 

We have made three images for each source, with increasing res- 
olution. To increase the sensitivity to partially extended sources, 
we have made naturally-weighted images from data with a cir- 
cular 30 % Gaussian taper at a (u, v) distance of 10M/1, yielding 
a resolution of 55.5x22. 1 mas 2 . For higher resolution and maxi- 
mum sensitivity, we have made naturally-weighted images using 
all data, yielding a resolution of 28.6x9.3 mas 2 . For maximum 
resolution we have made uniformly-weighted images using all 
data, resulting in a resolution of 20.8x5.9 mas 2 . All images were 
8192 2 pixels large, showing areas of 16.4x16.4 arcsec 2 for the 
tapered image and 8.2x8.2 arcsec 2 for the untapered images, re- 
spectively. Images were centred on the radio position derived 
from ATCA observations. 




19 20 21 



Fig. 2: Top left: Contour plot illustrating the relative orientation 
of the LCP (solid lines) and RCP (dashed lines) response of 
the VLBA antenna at Saint Croix. Contours are drawn at 95 %, 
90 %, 80 %, 50 % of the peak response. The series of crosses 
indicates the location of S462 in the antenna's primary beam 
during the observations, in intervals of approximately 6 min. Top 
right: The correction factor required to compensate for primary 
beam attenuation according to the location of the target, plot- 
ted as a function of observing time. The corrections for LCP 
and RCP can be quite different. Bottom left and right: The same 
plots for the Mauna Kea station. The series of crosses shows that 
the correction factors are markedly different from those at Saint 
Croix, because of the different parallactic angle. 



The naturally-weighted and uniformly-weighted images are 
still plagued by significant imaging artifacts, which we were un- 
able to remove. The sources appear to be irregularly extended 
predominantly in the north-south direction, in which the reso- 
lution is poor because of the low declination of the field. We 
attribute this to the fact that the visibilities which have the great- 
est resolution in the north-south direction are preferentially pro- 
duced near the rise and set of the source. At this time, the air- 
mass is high and these visibilities probably have comparatively 
large calibration errors. Another effect is that only relatively few 
long baselines are left after self-calibration on S503, because it 
is slightly resolved. These visibilities then make a disproportion- 
ately large contribution to uniformly-weighted images, reducing 
their quality. The loss of long baselines in the north-south direc- 
tion is reflected in a change of the ratio of the restoring beam 
axes, which increases from 2.51 (tapered images) to 3.1 (natural 
weighting) and to 3.5 (uniform weighting). 

3.4. Source detection and measurements 

The target source positions are relatively uncertain, in VLBI 
terms, as they were determined from arcsec-scale radio observa- 
tions that have a resolution 500 times lower than the VLBI obser- 
vations. To account for this, large VLBI images have been made 
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to increase the search space for source emission (see Sec. |3.3.3"| . 
We have used three criteria to automatically search for potential 
detections: 

1 . We have computed the SNR of the brightest pixel in the im- 
age (excluding a guard band of 100 pixels around the image 
edges which frequently have artificially high values arising 
from the imaging procedure). All images with SNR larger 
than 6 were treated as potential detections. 

2. The separation of the brightest pixels in the naturally 
weighted and tapered images was calculated. If this separa- 
tion was smaller than a beamwidth, a visual inspection was 
carried out. 

3. We have imaged the RCP and LCP data separately, using no 
taper and natural weighting. Since the two polarisations are 
processed by different electronics in the antennas, they pro- 
vide independent information about the sources. These im- 
ages were set to zero below a threshold, zl V2, with z — 5 <x 
and cr the noise in the Stokes I image, and then multi- 
plied. The resulting image is zero everywhere except at lo- 
cations where both the RCP and LCP images exceed zl V2. 
Compared to using a simple 5 <x cutoff in the Stokes I image 
this method is less susceptible to random noise peaks by a 
factor of more than 6 (see Appendix A). 

Sources matching at least one of these criteria were inspected 
by eye to confirm or reject the detection. 

To measure flux densities of the observed sources requires 
making a choice about resolution. Many sources are significantly 
resolved, and almost all sources can be reasonably imaged using 
heavily tapered or untapered data, resulting in a variety of resolu- 
tions. We decided to quote flux densities measured from images 
tapered at 10M/1, because at this resolution all sources (with the 
exception of S393) can be well approximated by a Gaussian. 
S393 has been fitted with two Gaussians, and the parameters 
listed in Table [3]represent their combined values. 

Flux densities were measured by fitting a Gaussian to the 
image, starting at the location of the brightest pixel. The most 
accurate results were obtained when the shape of the Gaussian 
was fixed, but its position angle kept as a free parameter. In very 
low-SNR cases such as ours, free fits frequently result in a gross 
overestimate of the flux densities. 

Flux density errors were estimated as follows. We have iden- 
tified three potential sources of error: (i) the a priori amplitudes 
in the VLBI observations are commonly assumed to be to be 
correct to within 10 %; (ii) the formal errors determined by the 
Gaussian fits; and (iii) the 10 % error estimate as a consequence 
of primary beam correction (which probably is conservative, in 
particular for sources near the pointing centre). These three com- 
ponents are added in quadrature to yield the total flux density 
error. 



Condon ( 1997 1 gives the equations for errors of Gaussian 



fits. In our case, the pixels are correlated, in fact, the area over 
which pixels are correlated is approximately one resolution el- 



ement in the images, or one beam area. For this case, Condon 
fi997\ gives 



fi 2 (l) ^(A) 



I 2 



A 2 



{ VM^m ) 



i 2 (9„ 



7 A/ 



2 



(3) 



where \i denotes a variance of a parameter returned by the 
fitting procedure, / the volume of the Gauss (the integrated flux 
density), A its amplitude (the peak flux density), and 9^, 9m, and 
9 m denote the FWHM of a circular Gaussian convolving function 



and the major and minor axes of the fit, respectively. In our case, 
~ 9m9„„ hence the equation simplifies to 



H\I) fi 2 (A) 



r- 



A 2 



fi z (9 M ) pL l {9 m ) 



9 2 



9 2 



(4) 



which we use as an estimate of the integrated flux density 
error arising from the fitting procedure. 

4. Results 

A summary of the image properties, source flux densities, and 
ancilliary data is listed in Table [3] 

4.1. Images 

Contour plots of the 20 detected sources and the one candidate 
detection are shown in Figure [3] The three images per source, 
made with different weighting and tapering, are shown in a row. 
The left panel shows a naturally-weighted image made with a 
10M/1 tapering (restoring beam size 55.5x22.1 mas 2 ), the centre 
panel the untapered, naturally-weighted image (restoring beam 
size 28.6x9.3 mas 2 ), and the right panel shows the uniformly- 
weighted image (restoring beam 20.8x5.9 mas 2 ). Positive con- 
tours start at three times the rms level of the images and increase 
by factors of y2. One negative contour is shown at three times 
the rms. 

5. Discussion 

5.1. Statistics 
5.1.1. Detections 

One of the goals of this project was to identify AGN in a blind 
survey of radio sources, i.e., without selecting objects other than 
by requiring a radio flux density in the ATCA survey. Given 
that the sensitivity of the VLBI observations varies by a factor 
of more than two over the observed area, the detection statis- 
tics are difficult to analyse. Also, our observations are biased 
towards compact objects, and this bias is a function of flux den- 
sity. Fainter sources will only be detected in VLBI observations 
if they are increasingly compact, whereas strong sources such as 
S393 (the second-brightest in our sample) can be detected if a 
fraction as low as 10 % is contained within a compact core. This 
means that the source types that our observations are sensitive to 
is a function of flux density. 

To provide some analysis of the detection statistics despite 
these shortcomings, we have divided our sample into three 
groups: sources with ATCA flux densities, S , of 10 mJy or more 
(bin 1, 7 objects), sources with 1 mJy < S < lOmJy (bin 2, 16 
objects), and sources with S < 1 mJy (bin 3, 73 objects). 

Out of the three bins, 5, 11 and 5 objects have been de- 
tected, respectively, corresponding to percentages of 71 %, 69 %, 
and 7 % (see Table [2). Most of the VLBA-detected sources (14 
sources) are compact, which is indicated by the ratio of VLBI- 
scale flux density to the ATCA-scale flux density being consis- 
tent with unity (s ee Table [3). W e note, however, that the ATCA 
flux densities by Norris et al. (20061 are probably affected by 
clean bias, which tends to underestimate the flux densities of 
weak sources. For a very similar observation, Middelberg et al. 
(2008a) show that clean bias can decrease the flux density by 
up to 5 % for sources with SNR smaller than 20. Hence the ra- 
tio of the VLBI-scale flux density to ATCA-scale flux density is 
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biased towards higher values, in particular for sources with an 
ATCA flux density of less than 1 mJy, and even exceeds unity in 
some cases. 



Table 2: The number of observed radio sources grouped into 
three flux density bins. A^atca is the number of radio sources 
found with the ATCA in a given flux density range, Nvlba is 
the number of sources detected out of these, %vlba is this num- 
ber expressed as a percentage of Najca, AGN is the number of 
known AGN prior to our observations, and nAGN is the number 
of new AGN we have identified. 



flux density 


^ATCA 


-^VLBA 


%VLBA 


AGN 


nAGN 


>10mJy 


7 


5 


71 % 


5 





1 mJy < 5 < lOmJy 


16 


11 


69% 


10 


5 


5 < 1 mJy 


73 


5 


7% 


19 


2 



5.1.2. New AGN 

It was noted in the introduction that a detection in a VLBI ob- 
servation is a strong indicator for AGN activity. To determine 
whether a source is already a known AGN we have used the in- 
formation from various sources listed in |Norris et al.| ( [2006| ) and 



the criteria by Szokoly et al. (2004i (using the X-ray hardness 



5.2.3 1 



ratio and X-ray luminosity, see Sec. 

Without the VLBI data presented here, the combined radio, 
IR, optical, and X-ray data were able to classify 36 sources 
(38 %) as AGN. There are 7 sources however, which have not 
been identified as AGN and are detected in the VLBI observa- 
tions. This number is 7 % of the total number of radio objects 
observed, but 35 % of the sources detected with the VLBA. 

All VLBI-detected bright sources in bin 1 were known to 
be AGN before, but out of the 16 medium-bright sources in bin 
2, only 10 (63 %) were previously known to be AGN, and our 
VLBA observations have added 5 (31 %). Among the 73 faint 
sources in bin 3, 19 sources were known AGN (26 %), and we 
have added 2 (4 %). These statistics suggest that VLBI observa- 
tions can make a significant contribution to the classification of 
galaxies below a flux density of 10 mJy, because brighter objects 
are already known AGN, and objects fainter than 1 mJy are only 
occasionally detected. This ability, of course, is a function of the 
sensitivity of the VLBI observations, which will soon improve 
significantly as the bandwidths of radio telescopes increase. 

Two sources, S423 and S443, have not been classified as 
AGN based on their optical and X-ray properties, but are de- 
tected in our VLBI observations. In these two cases the ancilliary 
data indicated a starburst or elliptical and a spiral galaxy, respec- 
tively. S423 is detected marginally in the uniformly-weighted 
image, with a resolution of 20.8x5.9 mas 2 . Its peak flux den- 
sity is 300 yt/Jy/beam, which indicates a brightness temperatures 
of 2.2 x 10 6 K. Hence it displays characteristics of a normal 
galaxy while at the same time exhibiting a brightness temper- 
ature clearly indicating an AGN. The VLBI detection of S443 
is offset from the centre of the galaxy and can therefore not be 
taken as evidence for an AGN; for further discussion on possible 
see Section 



5.1.3. Radio spectral index 

Radio spectral index can be used as an indicator for AGN activ- 
ity. When the spectral index is larger than -0.3 (S oc V) then 
the emission is deemed to come from self-absorbed synchrotron 
emission such as that found in the very compact regions of AGN 
and their jets. When the spectral index is smaller than -1.2 
then the emission is thought to come from extended structures 
formed by AGN in the high-redshift Universe. Between these 
values, a separation between AGN and star-forming activity can 



origins, 



5.3 



It has also been ignored in the follow- 
ing analyses which compare properties of VLBI-detected and 
VLBI-undetected sources. 



not be made. We have taken the data measured by Kellermann 
|et al. ( |2008[ ) and found that the median spectral index among the 
VLBI-detected sources is -0.7, whereas the undetected sources 
have a median of -0.8. In both cases, the scatter of the distribu- 
tion as characterised by the median average difference (MAD) is 
0.2. A Kolmogorov-Smirnov test (K-S-test) returns a p of 0.80 
that the two samples are drawn from the same parent population, 
indicating that there is no significant difference in spectral index 
between sources detected and undetected with VLBI. 

One might expect that the VLBI-detected sources show 
higher spectral indices than the undetected sources, because the 
VLBI observations require the emission to come from smaller 
regions, which then have a tendency to have smaller optical 
depths. We argue that since the linear resolution of our obser- 
vations is of the order of 280 pc (VLBI flux densities were mea- 
sured from images with a resolution of 55x22 mas and a typical 
redshift is 1) the emission from the observed regions is still dom- 
inated by optically thin synchrotron radiation. Unfortunately, the 
number of VLBI-detected objects is too small to test this hypoth- 
esis statistically. 

5.1.4. Redshifts 

We have searched the literature for redshift measurements of 
our sample. In total, we have found redshifts for 82 sources, of 
which 19 were detected in our VLBA observations. The median 
redshift of the detected sources is 0.98 with a MAD of 0.29, 
whereas the undetected sources have a median redshift of 0.67 
with a MAD of 0.41. Again the samples are not significantly 
different, with a p of 0.18 in a K-S-test. 

5.1.5. Infrared-Faint Radio Sources (IFRS) 

IFRS are mysterious objects which can be strong (tens of mJy) 
radio sources but do not have near-infrared counterparts de- 
tected by the Spitzer observatory as part of the Spitzer Wide-area 
InfraRed Extragalactic Survey (SWIRE). IFRS are unexpected 
because all known classes of galaxy at z < 2 are expected to 
be detectable given the IR survey sensitivities. See Huynh et al. 
( |20T0| , |Norris et aL] ( |20T0| ) and |Middelberg et al'| ( |2010| i and 
their references for recent work on these sources. Previous VLBI 
observations of six IFRS have resulted in two detections ( |Norris 
|et al.||2007| |Middelberg et al.||2008b| , supporting the idea that 
IFRS are AGN-driven. There are three radio sources in our sam- 
ple which have been classified as IFRS: S415, S446, and S506, 
which have catalogued arcsec-scale flux densities of 1.2 mJy, 
0.3 mJy, and 0.2 mJy, but none was detected. 



5.2. Co-located X-ray observations 

X-ray observations are claimed to be a very direct tracer of AGN 
activity (e.g., Mushotzky 2004). The radiation originates very 
close to the supermassive black hole, it can leave these regions 
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rather unabsorbed, in particular at energies of a few keV, and 
there is little contamination from other sources, such as stars. 

The surveyed region has been observed intensively with 
the Chandra X-ray observatory. In particular, a 940 ks expo- 
sure of a 0.11 deg 2 region as indicated in Figure [I] known as 
"the" Chandra Deep Fie ld South (CDFS), has r esulted in the de- 
tection of 346 sources (Giacconi et al. 2002} . The region has 
been reobserved by Luo et al. (2008 1, adding more than 1 Ms 



to the total exposure time and resulting in a central sensitivity 
of 7.1 x 10~ 17 ergcirT 2 s . The total number of X-ray sources 
found in this observation was 578. Furthermore, a 240 ks ex- 
posure of a 0.3 deg 2 region, known as the Extended Chandra 
Deep Field South (ECDFS), covers entirely the area observed 
here (Lehmer et al. 2005 ). This observation, with a maximum 
sensitivity of 3.5 X 10 _ltl ergcm _2 s _1 , yielded the detection of 
762 point sources. All observations were carried out in the en- 
ergy range 0.5 keV-8 keV and provide information about source 
flux densities in a soft band (0.5 keV-2 keV) and in a hard band 
(2keV-8keV). We have cross-matched our detections to these 
catalogues. 

5.2.1. Identifications in the CDFS 

The region indicated with the medium circle in Figure[T]contains 
23 radio sources, of which 6 were detected by our VLBI obser- 
vations. A search for counterparts of the VLBI-detected sources 
in the catalogues by Luo et al. (2008 ) returned 5 matches. There 
is only one VLBI-detected source, S474, which clearly lies in 
the area covered by the X-ray observations, but is not detected 
in the X-ray. Whilst it lies towards the edge of the area covered 
by the 2 Ms exposure, the integration times of the sources sur- 
rounding S474 suggest that its location has been observed for 
around 1.3 Ms. A further 4 VLBI-detected sources have X-ray 
counterparts outside the most sensitive region (see Figure ph. 



5.2.2. Identifications in the ECDFS 



Repeating the search in the catalogue by Lehm eret al.| ( |200"5] l, 
which covers an area that includes all sources observed by us, 
returned 12 sources, 8 of which were identical to sources dis- 
covered in the search of the CDFS data. Hence only one VLBI 
source was detected in the deep CDFS observations but not in the 
shallower ECDFS observations. In total 7 VLBI sources have no 
X-ray counterpart at all. 

This result sheds an interesting light on X-ray observations 
as a means of filtering AGN: at a sensitivity level correspond- 
ing to integration times of between 240 ks and 2 Ms the fraction 
of VLBI-detected AGN with X-ray counterparts increases from 
12/20 = 60% to 5/6 =83%. 

5.2.3. Hardness ratios 
The hardness ratio 



HR 



//cts - S cts 
//cts + S cts ' 



(5) 



where //cts and S cts are the net counts in the hard and soft 
band, can be used as an indicator for the degree of obscuration 
of an object. Soft X-ray spectra (with HR < -0.2) generally in- 
dicate that objects are unobscured and that, in the case of AGN, 
one is looking at a type 1 object. Harder spectra (HR > -0.2) 
indicate obscured lines of sight as would be caused by a dusty 
circumnuclear torus and can therefore be taken as indicators for 




-27°50' 



3 h .'ia :, ':!0* 33 m 00' 32 m 30' 32 m 00' 

Right Ascension (J2000) 

Fig. 4: Diagram illustrating the locations of X-ray sources in the 
CDFS (crosses), in the ECDFS (squares) and the sources de- 
tected in our VLBI observations (circles). 



type 2 objects. The classification based on hardness ratio alone 
can be improved when the total X-ray luminosity is also taken 

( |20Q4l >. 



into account, as has been demonstrated by|Szokoly et al. 



Sources with X-ray luminosities below 10" W (10^ erg s _1 ) are 
treated as non-AGN, sources with luminosities between 10 35 W 
and 10 37 W are considered as AGN and sources with X-ray lumi- 
nosities in excess of 10 37 W are classified as QSO. The hardness 
ratio is then used to further distinguish between type 1 and type 
2 objects. 

We have calculated the hardness ratio for all 13 VLBI- 



detected sources with X-ray counterparts. Where both Luo et al. 
( |2008| l and |Lehmer et al.| ( |20~0 5 ) had made a detection we have 
adopted the measurements of the former observation since it is 
more sensitive and its errors are smaller. Also note that in the 
case of sources with very few counts a rigorous error estimate es- 
sentially requires a re-analysis of the raw data ( |Park et aL 2006), 
which is beyond the scope of this paper. We therefore do not give 
an error for the hardness ratios. For all 14 objects redshifts were 
available in the literature, hence the total X-ray luminosity could 
be calculated and the Szoko ly~et al.| ( |2004) l criteria were used 
for classification. The scheme resulted in 3 non-AGNs (23 %), 2 
and 3 type 1/2 AGNs (15 % and 23 %), respectively, and 5 type 
1 QSOs (39 %, no type 2 QSOs). Out of the VLBI-undetected 
sources, 25 had sufficient data available for this classification, 
resulting in 10 non-AGN (40 %), 4 and 1 1 type 1/2 AGNs (16 % 
and 44 %), and no QSOs. See Figure|5]for an illustration of these 
statistics. It is difficult to read a trend off this distribution, how- 
ever, the lack of undetected type 1 QSOs is striking. Whilst only 
~ 10 % of optically selected type 1 QSOs are radio-loud, a much 
larger fraction of X-ray selected type 1 QSOs appear to be radio- 
loud. 

Part of the cause of this effect is that all 5 type 1 QSOs are 
brighter than 1 mJy, where the fraction of VLBI detections is 
generally high (69 %). 
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Fig. 5: X-ray classification statistics. Black bars represent 
sources with a VLBI detection and grey bars sources without. 
The lack of undetected type 1 QSOs is striking (see discussion 
in text). 



VLBI-detected 
VLBI-undetected 



1 10 100 

AT CA 1.4 GHz flux density / mJy 



Fig. 6: Full-band X-ray counts as a function of 1.4 GHz flux 
density. Neither the VLBI-detected (asterisks) nor the VLBI- 
undetected radio sources (pluses) show a correlation with X-ray 
counts. 



5.2.4. Potential causes of the X-ray non-detections 



Since X-ray observations have 
good tracer of AGN activity (e.g., 



been claimed to be a very 
Mushotzky|[2004l |Brandt~& 



Hasinger 2005), the lack of X-ray counterparts to almost 1/3 of 
the VLBI-detected sources deserves a closer look. 

First, it must be noted that the VLBI-detected source without 
X-ray counterpart in general have lower flux densities in other 
bands, too: only one of these objects, S329, does have a mea- 
sured optical magnitude, and only one other, S519, has a mea- 
sured 24 fim flux density. One object, S545, does have no optical 
or infrared flux density at all. 

Second, there appears to be no correlation between X-ray 
and radio flux densities (Figure IS). Hence some objects will un- 
avoidably have no X-ray counterpart, just as some sources with 
X-ray detections have no VLBI detection. 

Third, some sources could be Compton-thick and be so ob- 
scured that they remained undetected in the X-ray observations. 
This happens at column densities of Nu » 1.5 x 10 24 cm~ 2 
( Brandt & Hasinger 2005 ), which reduces the X-ray flux by two 
orders of magnitude, even at high (> lOkeV) energies. 

5.3. Notes on individual sources 

- S329 This source appears to be mildly resolved in our VLBI 
image, and is undetected at 0.5 keV - 8 keV by Lehme fet al.| 
(2005). Its redshift is 1.00 ( |Mainieri et al.|2008) . 

- S331 This source has been classified as AGN because of its 
more than tenfold radio excess over the radio-infrared re- 
lation. Its value of ^24 is -L96 which clearly puts it into 
the AGN regime. Its X-ray properties ( Lehmer et al.||2005~ 



|Szokoly et al.|2 004) suggest that it is a type 1 AGN. 
S380 This compact radio object was classified as AGN based 
on its q24 value of -0.21 ( |Norris et al.|2006] ). Its X-ray prop- 
erties qualify it as a galaxy because its X-ray lumino sity in 
the 0.5 keV -8keV band is only L x = 4.6 x 10 31 W (|Lehmer 



et al.||2005| |Szokoly et al.||2004| i. It has a spectral ind ex of 



0.6 pCellermann et al.||2008j and a redshift of 0.02 ( |Woif| 



et al.|2004 i, making it the closest object in our sample. 
S393 This is the only VLBI-detected source with a pro- 
nounced extension in the VLBI observations. It is an X-ray 



source at a redshift of 1.07 ([Zheng et a l. 2004), and it has 
been identified as AGN by |Norris et al.| (|2006|l bec ause of 
its X-ray hardness ratio. The Szokoly et al. ( 2004 1 criteria 
result in a type 2 AGN classification, and its radio spectral 
index is -1.2 ( |Kellermann et al.] [2008). Only a very small 
fraction (5 %) of its arcsec-scale flux density is recovered in 
the VLBI observations. At its redshift, the angular separation 
of 32 mas of the two components seen in our VLBI images 
translate to a projected distance of 263 pc, and is therefore 
in the regime of GPS sources (and is of a size comparable to 
the narrow-line region). 

S404 The VLBI images of this source display almost all 
the flux density found in the ATCA observations. It has a 
compact morphology and its X-ray properties classify it as a 
type 1 QSO ( |Luo et aLpOOB] |Szokoly et al.|2004) . It has a 
spectral index of 0.0 ( |Kellermann et al. 2008), indicating a 
compact, self-absorbed radio source, and its redshift is 0.54 



(Norris et al. 2006). The HST images show an unresolved 



object with diffraction spikes, indicative of a point-like op- 
tical source. Its q24, however, is 0.31 which is only a mild 
radio excess over the radio-infrared relation. 
S411 This is a compact source - all arcsec-scale emission 
is recovered in the VLBI image. Its hardness ratio of -0.44 



and X-ray luminosity qualify it as a type 1 QSO (Luo et al 
|2008|JSzokoly et al. [2004) 1 and it is located at a redshift ol 

has a spectral index of -0.5 



1.61 (Afonso et al. 



(Kellermann et al. 2008 ), and it is a featureless, round ob- 
ject in the HST ACS images ( |Giavalisco et al.|2004) >. 
S414 This is a compact source - the VLBI observations re- 
cover 95 % of the flux density measured with the ATCA 
( |Norris eP al. 2006). It has ^24 = -0.29 which makes this 
an AGN. S414 has a redshift of 1.57 ( |Mainieri et aLl|2008] l 
and is the brightest X-ray source in our sample, with Lx = 
1 .5 x 10 38 W (0.5 keV - 8 keV), qualifying it as a type 1 QSO 
( |Lehmer et al.|2005| [Szokoly et al.|2004) . 
S421 The radio data for this source are consistent with a 
point source since it is unresolved in our VLBI image. It 
is not an X-ray source, and it has a spectral index of -0.6 
(|Kellermann et al.|[2008) l and a redshift of 0.13 ( |Mainieri| 
|et al.|2008) . ~ 
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S423 This compact radio source is classified as a normal 
galaxy using the scheme by Szokoly et al.| (|2004| ). It has a 
spectral index of -0.2 ( Kellermann et al. 2008 ) and a redshift 
of 0.73 (Afonso et al. 2 006| l. Afonso et al. report spectro- 
scopic evidence for star-forming activity in this galaxy, po- 
tentially merging with a nearby (0.5") object. They identify 



five objects closer than 30" with similar redshifts. Mobasher 
et al. |( p004) l classify this object as an elliptical galaxy, us 



ing photometry in 17 bands. Hence there is no evidence ex- 
cept the VLBI detection that this object harbours an AGN. 
Given its redshift it has a 5 GHz radio luminosity of 7.5 x 
10 23 WHz which is about an order of magnitude too low 
to classify it as a radio-loud object, and is more typical of 
Seyferts. 

S437 This source has been classified as AGN based on <?24 
and the literature ( |Norris et al.|2006) . Only 4 % of its ATCA 
flux density were recovered, and the VLBI images show an 
unresolved source. Yet it is a powerful X-ray source exceed- 
ing a luminosity of 10 37 W (0.5 keV - 8keV) and it has a 



hardness ratio of -0.48, indicating a type 1 QSO (Luo et al. 
2008] |Szokoly et al.||2004|). The redshift of this object is 



0.73 ( jAfonso et al.||2006 i. It is one of four objects covered 
by the Great Observatories Origins Deep Survey (GOODS, 
Giavalisco et al.|2004 1 with the Hubble Space Telescope. It is 



an unresolved object in all four bands, and diffraction spikes 
indicate a strong point-like component. 
S443 This compact object is the faintest VLBI-detected radio 
source in our survey, with an arcsec-scale flux density of only 
0.4 mJy (Norris et al. 2006| >. It is associated with a spiral 
galaxy, with a spectral index of -0.5 ( Kellermann et al. 2008 ) 
and a redshift of 0.076 (Afonso et al.||2006[). Its ^24 value 



of 1.05 is broadly consistent with the value of 0.84 found 
by Appleton et al. ( 2004 ) for the radio-infrared relation, and 
also its X-ray properties suggest that it is a galaxy without 
AGN. Afonso et al. ( 2006| l report that the X-ray emission 
comes from north of the galaxy nucleus, with indications of 
star-forming activity. Mobasher et al. ( 2004 ) classify it as an 
Sbc object based on 18-band photometry. 
We consider this a marginal detection, partly because the 
VLBI position places the VLBI source at the very edge of 
the visible perimeter of the galaxy, which can be treated as 
evidence against a detection (see Figure 17). On the other 
hand, a radio transient cannot be ruled out (see |Lenc et aL] 
2008). We also note that the radio luminosity of this source, 
5 x 10 21 W Hz -1 , is similar to one of the brightest radio SNe, 
SN1986J. We therefore consider it possible that S443 is a 
radio supernova. 

S447 This source has been classified as AGN based on q24 
and also is an X-ray source. It has a spectroscopic redshift of 
1.96 ( Norris et al.|2006) , where 1 mas corresponds to 8.5 pc. 
It is unresolved in our images, hence most of the VLBI flux 
density (20 % of the arcsec-scale flux density) comes from a 
region smaller than ~100pc. It is the second most luminous 
X-ray source in our sam ple of 96 sources, w ith Lx = 1.3 x 
10 38 W (0.5 keV - 8 keV, |Lehmer et al.|2005| ). Nevertheless 
it has a soft spectrum with HR = -0.37 which qualifies it as 
a type 1 QSO ( |Szokoly et al.|2004| l. 

S472 This is a very compact source since 93 % of the ATCA 
flux density are seen in the VLBI image. Based on its <?24 = 
-0.31 this is an AGN. There is conflicting information about 
its redshift: |Zheng et al. (2004) give a photometric redshift 
of 1.22, but |Le Fevre et al. (2004) measured a spectroscopic 
redshift of 0.55, which we deem more reliable. Its X-ray 
spectrum is hard between 0.5 keV - 8 keV and indicates an 



absor bed type 2 AGN (Leh mer et aLl|2005| |Szokoly et ah] 
[2004l >. 

S474 This is a compact source which displays all of its 
ATCA flux density in a low-resolution VLBI image. Its red- 
shift is 0.98 ( Le Fevre et al.| [2004) and its spectral index is 
-1.4 ([Kellermann et al. 2008)>. It is a faint infrared source and 
is undetected in the X-ray. This is the only VLBI-detected ra- 
dio source which remained undetected by the 2 Ms Chandra 
exposure ( |Luo et al.|2008 1 



S482 This is a compact radio source. It is also an X-ray 
source (Lehmer et al. 2005) with hardness ratio -0.33, but 
its luminosity is only L x = 5.2 X 1 33 W (0.5 keV - 8ke V), 
which then qualifies it as a galaxy ( |Szokoly et al.|[2004) l. It 
was classified by Norris et al. ( 2006[ ) as an AGN based on its 
<?24 of -0.61. It has a spectr al index of -0.5 (|Kellermann et"aL 



2008 ) and redshift of 0. 15 dNorris et al.|2006| l. 
S500 This compact radio object has an absorbed X-ray spec- 
trum with hardness ratio 0.6 fLehmer et al.|2003]l , and it 
is classified as a type 2 AGN by the Szokoly et al. (2004) 
criteria. It has a spectral index of -0.9 ( Kellermann et al7| 
|2008] > and a redshift of 0.73 ( jMainieri et al.||2008) . Ba"sed 
on ^24 = -0.21, this is an AGN. 

S503 This is the brightest source in our VLBI observations, 
with a flux density of (21 .2+3.0) mJy (compared to an ATCA 
flux density of 19.3 mJy). Hence this source is very compact, 
but our uniformly-weighted image indicates that its core is 
resolved into two components with a separation of around 
lOmas (76 pc at the distance of S503; z = 0.81, IMainieri 



et al.||2008 1. It is an X-ray source indicating a type 1 AGN 



(Lehmer et al. 2005] |Szokoly et al.|2004) >. 
S517 This compact radio source is not a detected X-ray 



source, and it has a spectral index of -0.9 ( Kellermann et al. 
[2008] > and a redshift of 1.03 ( |Mainieri et al.|2008| >. TheeF 
tensions towards the south-east seen in the VLBI images are 
imaging artifacts and do not indicate a true extension. 

5519 This compact radio source has previously been clas- 
sified as AGN based on its <?24 = -0.42. It is not an X-ray 
source, and it has a redshift of 0.69 ( |Norris et al.|2006| ). 

5520 This source is mildly resolved. From its arcsec-scale 
flux density of 3.7 mJy (2.5 + 0.4) mJy could be recov- 
ered in our lowest-resolution image. There is a hint of a 
~25 mas extension towards the south-west in the naturally 
and uniformly-weighted image which at the redshift of the 
source of 1.4 (Wolf et al. 2004) corresponds to around 
200 pc. It has only a very faint IR counterpart in the co- 
located Spitzer image, and is not detected in the ECDFS 
survey by Lehmer et al. ( 2005 ) (it is outside the GOODS 
region). 

S545 This compact radio source also does not exhibit X-ray 
emission, and neither a spectral nor a redshift are available. 



6. Conclusions 

Wide-field VLBI observations of 96 sources in the Chandra 
Deep Field South have been carried out using a new correla- 
tion technique on a software correlator. The correlation was car- 
ried out only once, with high temporal and spectral resolution to 
avoid averaging effects, and the visibilities were phase-rotated to 
the 96 target fields inside the correlator, and then averaged. This 
approach is equivalent to a multi-pass correlation, but requires 
only a small computational overhead. It results in many small, 
normal-sized data sets, which can be calibrated using standard 
techniques. The results of this project are as follows. 
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Fig. 7: Left panel: HST B-band image of S443 (at the top) and 
S442 (bottom). Black contours (starting at 0.1 mJy and increas- 
ing by factors of V2) represent the ATCA 1.4 GHz flux density 
and the cross indicates the S443 VLB1 position. Right panel: An 
enlarged region of the HST image of S443 with enhanced con- 
trast and a cross indicating the VLBI position, which is at the 
very edge of the galaxy, but not outside it. 



- We have imaged all sources and have detected a fraction of 
21 %. The detection of one additional source is in the outer 
regions of the optical host galaxy, and could be a radio su- 
pernova. Most sources have flux densities of the same order 
as the arcsec-scale flux densities, and the radio-emitting re- 
gions therefore must be smaller than hundreds of pc (or less 
in some cases). Given their redshifts we can confidently in- 
terpret nearly all our VLBI detections as AGN. 

- A search of co-located, sensitive X-ray observations revealed 
that VLBI observations can identify AGN which have previ- 
ously been missed, even though X-ray observations are ex- 
pected to be a very good tracer of AGN activity. A total of 
7 sources were classified as AGN for the first time. Using 
X-ray data alone, only 10 out of the 13 detected sources with 
available X-ray counterparts were identified as AGN, miss- 
ing about one-quarter. 

- The VLBI core in the star-forming galaxy S443, which is 
one of the lowest redshift galaxies in our sample, has a radio 
luminosity of 5 x 10 21 WHz consistent with that of a su- 
pernova. The VLBI core in the galaxy S423 at a redshift of 
0.73 is surprising, because multiband photometry and X-ray 
data do not indicate the presence of an AGN. 

- Surprisingly, the VLBI detections include every X-ray de- 
tected type 1 QSO in our field, whereas only about 10% 
of optically-selected type 1 QSOs are radio-loud. Therefore, 
either X-ray observations preferentially select radio-loud 
QSOs, or even radio-quiet QSOs have compact cores de- 
tectable with VLBI. 



7. Appendix A: source finding using images of 
separate polarisations 

In large images, the probability of finding a random noise peak 
exceeding a given threshold can become significant. However, 
if the data are split in half and searched for emission indepen- 
dently, and if the positions are then compared between these two 
images, then the probability of finding a random noise peak can 
be lower than in the first case. We assume that the noise in the 
images has a Gaussian distribution, and that sources are unpo- 
larised - these assumptions are good approximations in our case. 

Suppose that the rms of a Stokes I =(RCP+LCP)/2 image 
is cr. Because the two orthogonal polarisations RCP and LCP 
are processed by independent antenna electronics the noise in 
images using only RCP or LCP is independent and larger than 
in Stokes I by a factor of V2. Hence locations in the Stokes I 
image which have an SNR of 5 have an SNR in the RCP and 
LCP image of5/V2 = 3.54. 

The number of random noise peaks, N', in an image above a 
threshold, z, is 

N' - N x p(z) (6) 
where 

P(z) = 1 - -L= P exp(4' 2 )* (7) 

and N is the number of resolution elements in the image. 
For example, for z = 5, p(z) evaluates to 2.87 x 10~ 7 , and in 
images which have 8192 2 pix and a resolution element consists 
of 209 pix, there are N' = 8192 2 pix/209pix x 2.87 x 10~ 7 = 
0.092 resolution elements above a 5 cr cutoff. 

To search for sources with the same absolute flux density 
in images made from RCP and LCP only, one has to lower the 
search threshold to zl V2, which results in a much greater num- 
ber of random noise peaks above the threshold. However, an ad- 
ditional constraint can be imposed on the search now, because 
one requires that peaks in the RCP and LCP images occur at the 
same location. 

Randomly selecting one resolution element out of N from 
RCP and LCP each results in l/N chance coincidences. 
Selecting one resolution element in RCP and N' from the LCP 
image (the number of resolution elements above the cr/ V2 
threshold) therefore results in N' x 1/7V chance coincidences. 
Repeating that N' times then results in N' 2 /N chance coinci- 
dences. Hence the number of random noise peaks occurring in- 
dependently at the same location in RCP and LCP when a detec- 
tion threshold is set at z = 5 cr in the Stokes I image is 

N" = N' 2 /N = (N p{zl V2)) 2 //V = N p{zl ^2? (8) 

Using the numbers from the example above, p(z/ V2) evalu- 
ates to 2.05x 10~ 4 , and N" = 8192 2 pix/209pixx(2.05xl0- 4 ) 2 = 
0.014. 

In our case the number of random noise peaks matching 
these criteria is 0.092/0.014 = 6.6 times lower than a simple 
5 cr cutoff. 
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Fig. 3: Contour plots of the detected sources. Three images per source are shown in a row. Left column: naturally-weighted images 
made with a 10 Mi tapering (restoring beam size 55.5x22.1 mas 2 ); middle column: untapered, naturally-weighted image (restoring 
beam size 28.6x9.3 mas 2 ); right column: uniformly-weighted image (restoring beam 20.8x5.9 mas 2 ). Positive contours start at three 
times the rms level of the images and increase by factors of V2. One negative contour is shown at three times the rms. 
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3 16.365 16.360 16.355 16.350 16.345 
RIGHT ASCENSION (J2000) 
Cont peak flux = 1 .0648 E-03 JY/BEAM 
Levs = 2.441 E-04 * (-1, 1, 1.414, 2, 2.828, 4, 5.657) 
CONT:S519 IPOL 1382.490 MHZ S519 low pb.lCLOOl .1 



03 33 10.210 10.205 10.200 10.195 10.190 
RIGHT ASCENSION (J2000) 
Cont peak flux = 9.21 80 E-03 JY/BEAM 
Levs = 2.496E-04* (-1, 1, 1.414, 2, 2.828,4, 5.657, 
8, 11.31, 16, 22.63, 32, 45.25, 64, 90.51) 
CONT:S517 IPOL 1382.490 MHZ S517 uni pb.lCLOOl .1 



- -27 47 24.80 - 



; 6.365 16.360 1 6.355 1 6.350 16.345 1 
RIGHT ASCENSION (J2000) 
Cont peak flux = 8. 7450E-04 JY/BEAM 
Levs = 2.148E-04 * (-1, 1, 1.414, 2, 2.828, 4, 5.657) 
CONT:S519 IPOL 1382.490 MHZ S519 lull pb.lCLOOl .1 



'50 1 6.745 1 6.740 16.735 1 6.730 16.725 
RIGHT ASCENSION (J2000) 
Cont peak flux = 9. 9583 E-04 JY/BEAM 
Levs = 2.578E-04*(-1, 1, 1.414, 2, 2.828, 4, 5.657) 
CONT:S520 IPOL 1382.490 MHZ S520 low pb.lCLOOl. 1 



16.365 16.360 16.355 16.350 16.345 

RIGHT ASCENSION (J2000) 
Cont peak flux z -5.5918E-04 JY/BEAM 
Levs = 2.596E-04" (-1, 1, 1.414, 2, 2.828, 4, 5.657) 
CONT:S519 IPOL 1382.490 MHZ S519 uni pb.lCLOOl .1 




50 16.745 1 6.740 1 6.735 1 6.730 1 6.725 
RIGHT ASCENSION (J2000) 
Contpeakflux= 7.801 7E-04 JY/BEAM 
Levs = 2.249E-04 ' (-1, 1, 1.414, 2, 2.828, 4, 5.657) 
CONT:S520 IPOL 1382.490 MHZ S520 full pb.lCLOOl .1 



-28 00 1 6.00 



50 16.745 16.740 16.735 16.730 16.725 
RIGHT ASCENSION (J2000) 
Cont peak flux z 6.3352E-04 JY/BEAM 
Levs = 2.777E-04- (-1, 1, 1.414, 2, 2.828, 4) 
CONT:S520 IPOL 1382.490 MHZ S520 uni pb.lCLOOl. 1 



-28 00 16.00 



16.775 16.770 16.765 16.760 16.755 
RIGHT ASCENSION (J2000) 

Cont peak flux = 1.9300E- 03 JY/BEAM 

Levs = 3.055E-04 "{-1,1,1 114, 2, 2.828, 4, 5.657, 



03 33 16.775 



16.765 16.760 16.755 
RIGHT ASCENSION (J2000) 
Cont peak flux = 1.1194E-03 JY/BEAM 
Levs = 2.661 E-04 -(-1.1,1 114, 2, 2.828, 4, 5.657, 



16.775 16.770 16.765 16.760 16.755 
RIGHT ASCENSION (J2000) 
Cont peak flux - 6.8897E-04 JY/BEAM 
Levs = 3.291 E-04 * (-1 , 1 , 1.41 4, 2, 2.828, 4, 5.657, 



Fig. 3: Continued 
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CONTiS545 IPOL 1382.490 MHZ S545 low pb.lCL001.1 



CONT: S545 IPOL 1332.490 MHZ S545 full pbjCLOOIj 



- -27 53 32.60 



03 33 33.445 33.440 33.435 33.430 33.425 
RIGHT ASCENSION (J2000) 
Cont peak flux = 1.1148E-03 JV/BEAM 
Levs = 3.131 E-04 * (-1, 1, 1.414, 2, 2.828, 4) 



CONT:S545 IPOL 1332.490 MHZ S545 uni pb.lCLOOl .1 



03 33 33.445 33.440 33.435 33.430 33.425 33.420 
RIGHT ASCENSION (J2000) 
Cont peak flux = 8.6938 E-04 JY/BEAM 
Levs = 2.737E-04 * (-1, 1, 1.414, 2, 2.828, 4, 5.657) 



03 33 33.445 33.440 33.435 33.430 33.425 
RIGHT ASCENSION (J2000) 
Cont peak flux = 7.81 05E-04 JY/BEAM 
Levs = 3.322E-04 -(-1.1,1 -114, 2, 2.828, 4) 



Fig. 3: Continued 



Enno Middelberg et al.: Wide-field VLBA Observations of the Chandra Deep Field South 



Enno Middelberg et aL: Wide-field VLBA Observations of the Chandra Deep Field South 
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